Stability of the hydrogen atom of classical electrodynamics 



O 

o 

(N 
o 

Q 



Oh: 



Jayme De Luc£0 
Universidade Federal de Sao Carlos, 

Departamento de Fisica 
Rodovia Washington Luis, km 235 
Caixa Postal 676, Sao Carlos, Sao Paulo 13565-905 
(Dated: February 2, 2008) 

We study the stability of the circular orbits of the hydrogen atom in a version of classical elec- 
trodynamics that was historically overlooked. We introduce the concept of resonant dissipation, i.e. 
a motion that radiates the center-of-mass energy while the interparticle distance performs bounded 
oscillations about a metastable orbit. The stability mechanism is established by the existence of 
a quartic resonant constant generated by the stiff eigenvalues of the linear stability problem. This 
constant bounds the particles together during the radiative recoil. The condition of resonant dissi- 
pation predicts angular momenta for the metastable orbits in reasonable agreement with the Bohr 
atom. The principal result is that the emission lines agree with the predictions of quantum elec- 
trodynamics (QED) with 1 percent average error even up to the 40*'' line. Our angular momenta 
depend logarithmically on the mass of the heavy body, such that the deuterium and the muonium 
atoms have essentially the same angular momenta, in agreement with QED. Analogously to QED, 
our stability analysis naturally uses the eigenvalues of an infinite-dimensional linear operator; the 
infinite-dimensionality is brought in by the delay of the electromagnetic interaction. 



O 
Oh' 



> 

(N 
(N 
O 

O 
o 



X 



According to classical electrodynamics, a quasi- 
circular orbit of the hydrogen atom radiates energy. We 
show that this complex Poincare- invariant two-body mo- 
tion can damp the center-of-mass energy only, while 
the interparticle distance performs bounded oscillations 
about some circular orbits, a phenomenon henceforth 
called resonant dissipation. This dynamics emits radia- 
tion of a sharp frequency while the center-of-mass energy 
decreases to compensate for the energy damping, during 
the life-time of the emission process, i.e. a collective ra- 
diative recoil. The stability mechanism is established by 
the existence of a resonant constant that keeps the in- 
terparticle distance bounded during the radiative recoil 
and deformation. This resonant normal form exists be- 
cause of a resonance involving the stiff modes of the lin- 
earized delay equations of motion, i.e., the linear normal 
modes with eigenvalues of a very large imaginary mag- 
nitude 0. The results of this simple stability analysis 
contain no adjustable parameter, satisfy Maxwell's equa- 
tions and predict the magnitudes of the Bohr atomj2j 
with accuracy and qualitative detail. We calculate every 
circular emission line of hydrogen with a 1 percent error. 
In this paper we also introduce an electromagnetic setting 
slightly different from the most popular form of electro- 
dynamics with self-interaction and retarded fields. This 
most popular form is henceforth called the dissipative- 
retarded setting (DR). Here we take the equations of 
motion derived from Fokker's Lagrangian of the action- 
at-a-distance electrodynamics Ja i3|, i^], with advanced 
and retarded interactions, and we add the Lorentz-Dirac 
self-interaction force 5] on each particle. This new elec- 
tromagnetic setting is henceforth called the dissipative- 
Fokker model (DF). Both the DR and the DF settings 
satisfy Maxwell's equations everywhere and define covari- 
ant dynamical systems for the electrodynamics of classi- 



cal point charges. We experimented with the stability 
analysis of the several possibilities within classical elec- 
trodynamics, and DF is the most transparent setting for 
the resonant dissipation phenomenon introduced here. 
Our results advocate for the interest of DF; we predict 
the Bohr atom with great precision and qualitative de- 
tail. It is popular to use the word radiation as equivalent 
to dissipation, because DR electrodynamics is absolutely 
dissipative. There is nevertheless a subtle difference from 
Poyting's theorem inside DF to Poyting's theorem in- 
side DR; the energy flux is not necessarily dissipative in 
DF (because we have advanced and retarded Lienard- 
Wierchert far-fields). This allows the atom to receive 
energy from other atoms of the universe via the effective 
half/advanced plus half/retarded far-fields, a thermaliza- 
tion mechanism that is absent in DR. 

To understand how our dynamics is a solution of the 
equations of motion of classical electrodynamics, one 
needs to go beyond the simplifications suggested by the 
Galilei-invariant Coulomb problem. In quantum elec- 
trodynamics (QED) the semiclassical Bohr orbits of hy- 
drogen correspond to excited quantum states that decay 
to the ground state in a life-time of about 10® turns[3|- 
The analog of this decay process is troublesome in DR 
and subtle in DF. Disregarding the singular coupling 
to the center-of-mass coordinate is a naive use of the 
equations of motion of DR, which predicts that the ra- 
diated energy comes out of the electrostatic potential 
energy. In this process the interparticle distance de- 
creases while the rotation frequency increases by a fac- 
tor of two during the life-time, making the emission of 
a sharp line impossible. It is very instructive to esti- 
mate this energy dissipation along hypothetical circu- 
lar orbits directly from the Lorentz-Dirac self-interaction 
forced- The complexity disregarded by this simplifica- 
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tion stems from three main sources: (i) The group the- 
oretical complexity of Poincare-invariant two-body dy- 
namics; the two no-interaction theorems^ limit Hamil- 
tonian and local Lagrangian descriptions of two-body dy- 
namics to straight-line free motions only. These negative 
theorems are sound warning that the Galilei-invariant 
Coulombian dynamics, with trivial and separable center- 
of-mass motion, is a bad approximation to relativistic 
two-body dynamics, (ii) The delay equations of mo- 
tion bring in an infinite-dimensional dynamical system 
and needs an initial function as the initial condition; 
addition of the non-runaway condition leads to a well- 
posed mathematical problem 0. The stability analysis 
of such dynamics reveals stiff eigenvalues with an arbi- 
trarily large imaginary part, as long as integer multiples 
of TT (an arbitrary integer number appears naturally) T| . 
(iii) The generic solution to a stiff differential equation is 
not slowly varying everywhere and typically it has jum^; 
see for example the equations of Lienard type in Ref . s\ . 
Another example is the small parameter multiplying the 
highest derivative of the Lorentz-Dirac self-interaction 
forceQ. After Dirac's derivation of a finite force of self- 
interaction within classical electrodynamics 5] , Eliezer 
considered a model two-body motion with Coulomb in- 
teraction plus self-interaction'^. A simple physical con- 
sequence transpired, that with self-interaction the elec- 
tron cannot fall into the proton even along unidimen- 
sional collision orbits. This surprising result is a warning 
that self-interaction is a singular perturbation over the 
Coulombian dynamics and also suggests that colinear or- 
bits are the natural attractors of the dissipative dynamics 
(a ground state with zero angular momentum). Along 
such orbits the heavy particle (the proton) moves in a 
non-Coulombian way and the self-interaction provides a 
repulsive mechanism that avoids a collision. Disregard of 
the stiff nature of the equations of classical electrodynam- 
ics has produced several conundrums, and the present 
work is an attempt to recognize the correct qualitative 
behavior of the solutions of classical electrodynamics. 

In Ref. |0| we studied the classical helium atom in the 
DR setting. Because the very special circular orbits of he- 
lium do not radiate in dipole inside DR, the physics of the 
emission of a sharp line is easy to understand even in the 
DR setting^3|- There is already an interesting resonant 
structure for these circular obits of helium within the 
Darwin approximation (a low-velocity approximation to 
Fokker's Lagrangian) [llJ,|ll| The predictions of the Dar- 
win approximation were tested numerically in Ref. |l2l |. 
finding non-ionizing states in the atomic magnitude, as 
predicted in Ref. [lOj . The generalization to the circular 
orbits of the hydrogen atom in the DR setting is trou- 
blesome; circular orbits of hydrogen are stable within the 
Darwin approximation and radiate in dipole, such that 
a sharp emission line seems like a conundrum for classi- 
cal electrodynamics in the DR setting. In this work we 
come to grips with the physics of the hydrogen atom by 



changing to DF and by understanding that the correct 
description comes only after the inclusion of the nuclear 
motion (i.e. a two-body system with delay and dissi- 
pation). The resonant dissipation radiates mostly the 
center-of-mass energy while a slow deformation of circu- 
larity and a collective recoil take place. After a possibly 
long life-time, this resonant dissipation produces enough 
deformation to destabilize the orbit, i.e. only metastabil- 
ity is possible (this is also what QED says about circular 
orbits). 

Our perturbation theory starts from the unperturbed 
circular orbit of the isolated electromagnetic two-body 
problem derived from Fokker's Lagrangian, as explained 
in Refs.P, |^ Q. A motion of the particles in concen- 
tric circles with a constant angular speed and along a 
diameter is a solution of these equations of motion. This 
is because the symmetric contributions from future and 
past generate a resulting force that is normal to the ve- 
locity of each particle. We henceforth call this orbit the 
circular orbit [3, H, 0. The quantities for this uniform 
circular motion are given in Our electron has mass 
TOi and travels a circle of radius ri while the proton has 
mass m2 and travels a circle of radius r2 ■ We henceforth 
use a unit system where the electronic charge is e = — f 
and the speed of light is c — 1 . The circular orbit is de- 
fined by the retardation angle 6 that one particle turns 
while the light emanating from the other particle reaches 
it (one light-cone distance away). The frequency of an- 
gular rotation is given, to lowest order in 9, by Kepler's 
law 

n = ne^ + ... (I) 

where fj. = mi 7712/(7711 + 7^2) is the reduced mass. For 
shallow-energy atomic orbits, the angular momentum 1^ 
is given to first order by = which is of the 

order of the inverse of the fine structure constant, = 
137.036. To simplify the stability analysis, we transform 
the particle coordinates from Cartesian to gyroscopic co- 
ordinates 

Xk+iyk = exp(7rit)[dfc -f 7/4,], (2) 
Xk-iyk = exp(-7f7i)[dfc -I- ^fe], 

where rjf. and are complex numbers defining the per- 
turbation of the circularity. The coordinates of the elec- 
tron are defined by Eq. Q with k = 1 and di = ri while 
the coordinates of the proton are defined by Eq. (0) 
with k — 2 and c?2 = —r2 (at the same time in the iner- 
tial frame the particles are in diametrically opposed posi- 
tions along the unperturbed orbit Q)- We also need here 
the stability along the perpendicular z direction, a linear 
problem that is uncoupled from the xy stability and eas- 
ier to work out. The linear stability analysis proceeds like 
in Ref. Q], by substituting the circular orbit plus a per- 
turbation, Eq. into Fokker's Lagrangian, expanding 
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this Lagrangian to quadratic order and taking the Euler- 
Lagrange equations to Unear order 0. Because we have 
a delay system at hand, the normal mode condition (the 
characteristic equation) has infinitely many roots. This 
characteristic equation involves hyperbolic functions that 
can become arbitrarily large and are insensitive to an ar- 
bitrarily large imaginary part of the eigenvalue [J. A 
standard technique of delay analysis is to keep only the 
largest powers of the eigenvalue plus the hyperbolic terms 
in the characteristic eauation|l4l|. This large-imaginary- 
part stiff limit is performed in Ref. 0] for the equal-mass 
two-body system (see Eq. (15) of Ref. 0). Here we gen- 
eralize this large-imaginary-part limit to the arbitrary- 
mass case and we computed the characteristic equation 
with a symbolic algebra software. We define the normal 
mode eigenvalue by i.e. every coordinate pertur- 

bation oscillates in time as exp(XQt / 9) ( A is an arbitrary 
complex number). The limiting form of the characteristic 
equation for the isolated different-mass case is 



(^)cosh2(A) = l, 



(3) 



where fj, is the reduced mass and M = mi+m2. Both the 
planar and the perpendicular linearized equations have 
the same limiting characteristic Eq. ^ along circular 
orbits. For hydrogen (fi/M) is a small factor of about 
(1/1824). It is important to understand the structure of 
the roots of Eq. Q in the complex A plane, specially 
for 9 of the order of the fine structure const^ant. It is 
easy to see that the very small parameter ~ 10"^'^ 
multiplying the hyperbolic cosine on the left-hand side 
of Eq.Q determines that a = \ Re(A)| ~ \n{,J^). For 

the first 13 excited states of hydrogen this a is in the 
interval 14.2 < \a\ < 18.2. The imaginary part of A can 
be an arbitrarily large multiple of tt, such that the general 
solution to Eq. Q is 



X — ±{a + inq), 



(4) 



where q is an arbitrary integer. The plus or minus sign 
of Eq. is related to the time-reversibility of the iso- 
lated two-body system, a symmetry that is broken by 
radiation. 

Next we include the dissipation of the DF electrody- 
namics, i.e. the Lorentz-Dirac self-interaction, a calcula- 
tion performed by adding the self-interaction force to the 
Lagrangian equations of motion of the isolated system. 
Here we give only the limiting form of the characteristic 
equation 



1 ./xr 



3^^A+-(^)sinh(2A) + . 



(l-^)(^)cosh2(A) = l- 

(5) 

where g = —7 for the planar stability and g = 1 for 
the perpendicular stability. We give the terms up to 



0(1/A^) because these are essential for the mathematical 
phenomenon presented. The coefficient g of the term of 
order 1/A^ is the only difference between the planar and 
the perpendicular characteristic equations up to 0(1/A^) . 
Further terms of size 1/A^ and 9*X^ produce only small 
quantitative changes for a in the atomic range and will 
be given elsewhere. The linear term on the right-hand 
side of Eq. © with the 2/3 coefficient is due to the 
self-interaction force. This dissipative term breaks the 
time-reversal symmetry of Eq. and the roots of Eq. 

no longer come in plus or minus pairs. In the follow- 
ing we implement the necessary condition for a convex 
analytic function of the interparticle distance to exist. 

In the same way already used in Ref. [lo| , this condition 
should involve the perpendicular direction z as well. This 
is because if the atom is to recoil like a rigid body while 
it radiates, one expects the orbital plane to nutate and 
precess, typical gyroscopic motions of a rigid body. To 
introduce the main idea simply, we assume that the non- 
runaway condition restricts the physics to a finite number 
of linear modes. We take a perpendicular normal mode 
and a planar normal mode of Eq. Q , of eigenvalues A^ 
and Xxy respectively. Suppose we can find these modes 
such that Xz + Xxy + A* -I- A* = 0, which is the necessary 
condition for a quartic normal form. i.e. 



Re(A^ + Xxy) = 0. 



(6) 



The coordinate of the planar normal mode is a linear 
combination of the four 77^ gyroscopic coordinates: u = 
o-ikVk +^ifeCfe , while the coordinate of the perpendicular 
z normal mode is Z = bizi + b2Z2- Because Fokker's 
Lagrangian is real. A* and A*j, are also solutions to Eq. 
(O with complex conjugate normal mode coordinates. 
Using these normal mode coordinates and Eq. lO, one 
can show that the following quartic form is a constant of 
motion up to higher order terms [llj: 



C : 



\Z\ 



(7) 



Notice from Eq. (0) that |up is a quadratic function of 
the planar coordinates, such that the resonant constant 
of Eq. (O limits the excursion along the z direction times 
the planar distance to the circular orbit, essentially lock- 
ing the dynamics to the circular orbit. This necessary 
condition and the continuation of the leading term |(7|) to 
an asymptotic series is discussed in Ref. llj. The les- 
son of the above construction is that if the interparticle 
distance is to be bound by a resonant constant while the 
dissipation goes on, then we must be able to find a pair 
of roots to Eq. such that Eq. © holds. This implies 
that Eq. lO has a pair of almost-symmetric roots of the 
form 



A 



xy 

A. 



(cr + nqi + zei), 
-((T + nqi + 162), 



(8) 
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where ei and €2 can be assumed real if Xz + Xxy is purely 
imaginary by absorbing the real part of the perturbation 
in the definition of a. The above-defined root-searching 
problem is well posed and for each integer q conditions 
(O and (jS)) determine 6* as a function of q , i.e., 9 must 
be quantized! An asymptotic solution to condition (jS)) 
can be obtained by expanding Eqs. lO up to quadratic 
order in ei and €2 while treating a as an approximate 
constant. This approximation determines the following 
discrete values for 



6(7r2g2 



a{'K'^q'^ 



and 



7r2g2)2 



(9) 



(10) 



According to QED, circular Bohr orbits have maximal 
angular momentum and a radiative selection rule ( AZ — 
±h) restricts the decay from level fc + 1 to level k only, 
i.e. circular orbits emit the first line of each spectro- 
scopic series (Lyman, Balmer, Ritz-Paschen, Brackett, 
etc.), henceforth called the QED circular line. We have 
solved Eqs. (jHJ and (O with a Newton method in the 
complex A plane. Every angular momentum Iz — 0^^ 
determined by Eq. ® has a value in the correct atomic 
magnitude, but the subset of Table 1 has frequencies 
wdf surprisingly close to the QED lines. These lines 
are when the integer q is approximately equal to an in- 
teger multiple of the integer part of 2a-. The reason for 
such selection rule is that from the resonances satisfy- 
ing the necessary condition, only some have jup depend- 
ing on the translation-invariant quantities — ^2) a-nd 
(rji — 772) to allow a recoiling translation 16J. In our de- 
scription the emission mechanism is at a frequency equal 
to the orbital frequency f2 corrected by the frequency 
of the complex amplitude uZ defined above Eq. |(7J), 
as we explain below. The numerically calculated angu- 
lar momenta Iz = for this select subset are given 
in Table 1, along with the orbital frequency in atomic 
units {I37^n)/^i = (13761)3, the QED first frequency of 
the series in atomic units wqed = 5(1^ ~ (FFTp")' ^^'^ 
the frequency predicted by the dissipative Fokker model 
WDF = (1376')3 -t- m^d'^iei - 62). We list only the first 
13 lines, which are the experimentally observable, but 
we tested the agreement of the numerical calculations of 
the Newton method with up to the 40*'* circular line pre- 
dicted by QED. Beyond that, the asymptotic formula @ 
shows that the agreement is essentially for any integer k 
because substitution oi q — [2a]k into Eq. ^ yields 



-1 = J^^3/2^^ 



137.9fc, 



(11) 



to be compared with the 137.036 of QED. The agreement 
for any integer k suggests that Eq. is equivalent to 



Schroedinger's equation (both are eigenvalue problems 
for an infinite-dimensional linear operator, and linear op- 
erators with the same spectrum are equivalent) . 



Iz = 


(13761)3 


Wqed 


Wdf 


161.21 


6.137x10-1 


3.750x10-1 


3.688x10-1 


282.88 


1.136x10^1 


6.944x10-2 


6.806x10-2 


397.40 


4.097x10^2 


2.430x10-2 


2.470x10-2 


519.59 


1.833x10-2 


1.125x10-2 


1.113x10-2 


637.79 


9.911x10-3 


6.111x10-3 


6.054x10-3 


751.48 


6.059x10-3 


3.685x10-3 


3.719x10-3 


871.84 


3.880x10-3 


2.391x10-3 


2.392x10-3 


987.25 


2.672x10-3 


1.640x10-3 


1.654x10-3 


1109.18 


1.188x10-3 


1.173x10-3 


1.170x10-3 


1225.92 


1.139x10-3 


8.678x10-"^ 


8.694x10-4 


1343.20 


1.061x10-3 


6.600x10-1 


6.628x10-1 


1466.89 


8.146x10-1 


5.136x10-"^ 


5.102x10-* 


1585.51 


6.455x10-"! 


4.076x10"'! 


4.052x10-1 



Table 1: Numerically calculated angular momenta 
Iz — in units of e2/c, orbital frequencies in atomic 
units (1370)3, circular lines of QED in atomic units 
Wqed = ^Ip" — (k+i)^ ) emission frequencies 

of the DF model in atomic units wdf = (1376')3 + 

Last we calculate the frequency of the emission line. 
The treatment is a consequence of the above defined 
DF version of electrodynamics and differs from the 
usual electromagnetic model (DR) in non-trivial ways. 
We have consistently used the relativistic action-at-a- 
distance electrodynamics as a many-body electromag- 
netic theory, as first suggested in Ref. 13] . This no-field 
theory is based on a parametrization-invariant action in- 
volving two-body interactions only, without the media- 
tion by fields. The isolated electromagnetic two-body 
problem, away from the other charges of the universe, is 
a time- reversible dynamical system defined by Fokker 's 
action 



Sf = — J midsi — J m2ds2 (12) 
+ J J d{\\xi - X2\\'^)xi ■ X2dsidS2, 

where Xi{si) represents the four-position of particle i = 
1 , 2 parametrized by its arc- length Si , double bars stand 
for the four-vector modulus \\xi — X2\\'^ = (a^i — X2) ■ 
{xi — 2:2), and the dot indicates the Minkowski scalar 
product of four- vectors with the metric tensor g^^ (goo = 
li5ii = 522 = <?33 = —1)- The Euler-Lagrange equations 
of Lagrangian ifT^ were already used to derive our Eq. 
©. Surprisingly, the solution to the unusual equations 
of motion of Lagrangian (|12|l are still determined by ini- 
tial position and velocity, as proved in Ref.(j2j) for the 
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non-runaway solutions. Therefore we are dealing with a 
perfectly well-posed and causal system dressed in an un- 
familiar form, and we invite the reader to consult Ref . ■ 
In the many-body action-at-a-distance electrodynamics 
one needs a prescription to describe the incoherent in- 
teraction with the universe, because including every sin- 
gle particle of the universe in a pairwise many-body La- 
grangian like Eq. H12(l is impractical (even though it is 
precisely what we want to describe physically) . This pre- 
scription for the dissipative interaction with the universe 
must be a covariant prescription that satisfies Maxwell's 
equations. For example, the Wheeler-Feynman prescrip- 
tion of Ref. ^3 leads to the usual electromagnetic DR 
model. The prescription of the DP model is to add 
the Lorentz-Dirac self-interaction directly to the Euler- 
Lagrange equations of action ifT^ . 

In the DP model the interaction with a distant par- 
ticle involves half/retarded plus half/advanced Lienard- 
Wierchert interactions, as derived from a pairwise La- 
grangian of type (|12|l . We shall show that this does not 
preclude radiation of energy. Once DP satisfies Maxwell's 
equations everywhere, Foynting's theorem is valid for 
any true orbit that solves the stiff equations of motion. 
Thinking of Poynting's theorem as detached from the 
complexity of the equations of motion might suggest a use 
of this theorem to calculate the radiation of hypothetical 
orbits that are not even solution to the equations of mo- 
tion. For example, one can show that with half/retarded 
plus half/advanced fields a perfectly circular orbit does 
not radiate energy. But again, it is just as easy to show 
that a perfectly circular orbit is not a solution of the 
equations of motion with the self-interaction terms. Stiff- 
ness can prevent the dynamics from staying even for a sin- 
gle turn in the neighborhood of an orbit that would be a 
solution to the equation without the stiff terms |1] . A fa- 
miliar use of the retarded-only radiation fields of DR is to 
calculate the incoherent dissipation of a circular current 
by applying Poynting's theorem to the radiation fields 
of a hypothetical circular orbit. It is a trivial exercise 
to show that this produces the same result of integrat- 
ing the Lorentz-Dirac self-interaction times the velocity 
along that hypothetical circular orbit. Even though this 
approximation is fine for macroscopic currents restricted 
by a wire and forced by a battery, the stiff delay equa- 
tions of electrodynamics can jump wildly given unstable 
conditions. Therefore, one should see first if the solu- 
tion exists and if it is stable in the model and only then 
calculate the radiation for that precise stable orbit. The 
perturbations to circularity defined by Eq. (jSJ are nec- 
essary to stabilize the orbit, and as we show below, are 
essential as well to determine the radiated energy. The 
self-interaction force is a very small singular dissipative 
term along the unperturbed circular orbit, and the dissi- 
pation calculated over the perturbed orbit is completely 
changed by the gyroscopic perturbations, even if the final 
orbit stays near the circular orbit. After proper account 



of these perturbations, our metastable orbit radiates even 
with the half-retarded plus half-advanced fields. 

The half/retarded plus half/advanced radiation mag- 
netic field of the electron in the DP model (the far- 
magnetic field) is[l5| 



Brad — 



(a_ X n_) 



(a+ X n+) 



2(1 — h- ■ V-)^r 2(1 — n+ • i>_|_)^r ' 



(13) 



where v and a are the electronic velocity and accelera- 
tion, n is a unit vector from the electron to the observa- 
tion point, underscore minus indicates evaluation on the 
retarded light-cone and underscore plus indicates evalu- 
ation on the advanced light-cone. These two light-cones 
are defined hy t± — t ± (r — h± ■ x) where x is the vec- 
tor position of the electron measured from the center of 
the circular orbit. Once (i+ — t„) ~ 2r, whenever 2r is 
an integer multiple of the distance travelled by light in 
a circular period (one light-period), the first term of ex- 
pansion (|13|l along the precise circular orbit vanishes and 
the next term is the important one; a quadratic function 
of the orbital quantities. The same cancellation happens 
for the far-electric field whenever 2r is an odd-integer 
multiple of half the light-period, but for us here the far- 
magnetic field is enough. The next term of the magnetic 
expansion about the circular orbit is 



(1) 2(n • v){a X h) + {h ■ x){d x h) 

^rad = I • 



(14) 



We can estimate H14f) by noticing that along the n± — x 
direction of the unperturbed plane this quadratic term 
contains a product of the z perturbed coordinate times 
the X perturbed coordinate. Those perturbations are de- 
scribed by the u and Z nonrunaway normal modes al- 
ready used and explained above Eq. (0. Using the nor- 
mal mode conditions 9u — ClX^yU and 6Z = QX^Z to- 
gether with Eq. (|SJ| one can show that the quadratic form 
uZ is a complex amplitude that oscillates harmonically 
with the beat frequency (Xxy + Xz)fl/0 = (ei — 62)^1/6'. 
This complex amplitude is precisely the square-root of 
our resonant-normal form (jjj), which is a constant be- 
cause it is the modulus of the complex number uZ = 
VC exp[i(ei — e2)Ot/0]. Translating the u mode to Carte- 
sian coordinates with Eq. we obtain the following 
approximation to i?, 



(1) 



ni 2uZ , „ , 

Kld^—^Mm. 



(15) 



and therefore the frequency of the emission line is equal 
to SI plus the frequency of uZ 



wdf = si + (ei - e2)ri/0. 



(16) 



with S7 given by Eq. Notice that the emitted fre- 

quency of the DP model is naturally different from the 
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orbital frequency. Among all possible settings of electro- 
dynamics, only DF has this possibility. The fact that 
the emission frequency of hydrogen is different from the 
orbital frequency is another famous conundrum for the 
usual DR electromagnetic setting. The emission fre- 
quency of Eq. 1)161 ) contains differences of eigenvalues of 
the linear infinite-dimensional operator ((SJ and is strik- 
ingly similar to the Rydberg-Ritz combinatorial principle 
of quantum mechanics for the emission lines. 

To discuss the width of the emission lines, we notice 
that for small amplitudes the resonant normal form con- 
strains the deformation while allowing a radiative recoil. 
Even though the first term of the resonant normal form 
is a function of the relative separation only, a higher- 
order term might not be so, and at some point the recoil 
breaks the metastable orbit down. The construction of 
this next term shall be left for future work. In the dy- 
namical process of resonant dissipation, the sharp line 
is emitted while the dynamics is locked in the neighbor- 
hood of the original orbit. We conjecture that when the 
metastable orbit breaks down, it can fall into the next 
metastable attractor; another circular orbit, or into the 
ground state. The stiff modes used in Eq. ((SJ describe 
a fast (stiff) time-scale of a frequency of the order of 
a/9 ~ 1400 times the orbital frequency, the time for a 
stiff jump of the dynamics. After this fast timescale the 
resonance essentially locks the dynamics to the neighbor- 
hood of the metastable circular orbit. 

The present simple approximation does not yet com- 
pete with the precision of QED, but it is surprising that 
recognizing the correct qualitative dynamical behavior of 
the equations of motion of classical electrodynamics can 
take us so far (the small corrections to Eq. jSJ will be 
given elsewhere) . The results are also beyond the Darwin 
approximation to Fokker's Lagrangian because the log- 
arithm is not analytic at 6* = 0. We see then that the 
physics of resonant dissipation is qualitatively beyond 
the non-ionization criterion previously used by us|l2l|. 
even though it contains it. Our stability condition forces 
the particles to stay bound during dissipation, a bound 
state condition analogous to the normalizability of the 
wave functions of Schroedinger's equation. Our stability 
analysis naturally involved an infinite-dimensional linear 
eigenvalue problem, Eq. ((SJ, which is essentially a PDE, 
like Schroedinger's equation. We conjecture here that it 
should be possible to construct an average linear opera- 
tor from Eq. ^ with g = 1 and Eq. Q with g = —7 
such that n + {ei — €2)^/9 is a difference between its con- 
secutive eigenvalues. That would be the Schroedinger's 
eigenvalue problem. Last, the angular momenta of Eq. 
()ll(l depend only logarithmically on the mass ratio times 



9 , such that the deuterium and the muonium atoms have 
essentially the same values for 9 (i.e., the same fine struc- 
ture constant, in agreement with QED). Notice that if 
the proton mass is set infinite, the quantized angular mo- 
menta become infinite as well, such that this logarithmic 
CT is a genuine two-body effect. 

Recognizing the correct qualitative dynamical behav- 
ior is an intuition that goes a long way here; it suggested 
we should explore the infinite dimensionality of the delay 
with stability analysis and make use of the stiff modes. 
These stiff modes provide a natural integer number to 
label the metastable orbits. The natural appearance of 
the infinite-dimensional operator (Eq. lO), and the large 
body of qualitative and quantitative agreement suggests 
that simple stability analysis inside the DF electrody- 
namics is a consistent new version of QED; one with an 
h given by Eq. Hll() . More work should be done to un- 
derstand the relation between those two theories. 

We thank Savio B. Rodrigues for discussions and for 
checking the numerical work in the complex A plane with 
MATLAB. 
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